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A  method  is  proposed  for  treating  linear  longitudinal  perturbations  in  one- 
dimensional  collisionless  plasma  diodes  with  a  uniform  plasma  region  and 
thin  electrode  sheaths.  The  method  is  comprehensive  in  that  it  allows  for 
very  general  equilibrium,  initial,  boundary,  and  external-  circuit  conditions,  c 
Upon  Laplace-transforming  the  Vlasov  and  Poisson  equations  in  both  space 
and  time,  appropriate  evaluation  of  all  pertinent  relations  leads  to  a  set  of 
2+2 n»  (n*  is  the  number  of  particle  species)  coupled  integral  equations  in  x 
and  v  for  the  following  quantities  (which  are  the  time  Laplace  transforms  of 
the  respective  physical  perturbations):  y,(w)  (external-circuit  current  den¬ 
sity),  E(z,u)  (electrostatic  field),  /*(»  >  0,w),  and  /*(»  <  0,w)  (velocity 
distribution  functions  of  the  plasma-bound  particles  at  the  left-  and  right- 
hand  plasma  boundaries,  reap.),  where  <r  is  the  species  index.  The  formal 
solution  of  these  integral  equations  and  the  inverse  Laplace  transformation 
(w  —  t)  are  discussed  in  general  terms.  In  particular,  it  is  shown  that  the 
intrinsic  eigenfrequencies  are  given  by  the  zeros  of  the  coefficient  determi¬ 
nant  of  the  integral  equations.  A  comparison  with  previous  treatments  is 
given,  and  it  is  concluded  that  extensions  of  the  method  proposed  to  more 


-f- 


general  systems  should  be  feasible. 
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L  INTRODUCTION  AND  SUMMARY 


In  the  first  part  (§1)  of  his  famous  1946  paper,1  Landau  derived  the  gen¬ 
eral  solution  to  the  initial-value  problem  for  small-amplitude  longitudinal  elec¬ 
tron  oscillations  in  an  infinite,  uniform  collisionless  plasma.  Any  perturbation 
o(z,  v,  t)  was  represented  as  a  spatial  Fourier  integral,  with  every  Fourier  com¬ 
ponent 

a(q,v,t)e~<1M  {q  real) 


evolving  in  time  independently  of  the  others.  (Here  and  henceforth  we  use  a 
notation  conforming  with  the  rest  of  this  paper.)  The  key  result  was  Landau's1 
Eq.  (10),  which  provided  an  explicit  expression  for  the  time  Laplace  transform 
of  the  potential-perturbation  amplitude  associated  with  a  given  wavevector  q. 


In  a  straightforward  extension  to  multi-species  plasmas,  Landau's1  equa¬ 
tion  (10)  can  be  re-written  in  the  form 

»# 


“>  “  7dm  5  /-»  *+“  h  •'  ” 


(1) 


where  q  is  real,  w  is  complex  (cf.  Subsec.  n.B),  V(q,u)  is  the  Fourier-Laplace 
transform  of  the  electrostatic-potential  perturbation  V(x,t),  a  =*  1 ,  ...,n„  is 
the  particle-species  index,  e*  is  the  electric  charge  of  a  particle  of  species  cr, 
f*{q,v)  is  the  Fourier  transform  of  ~f*(x,v)  (the  initial  perturbation  of  the 
velocity  distribution  function  of  species  <t),  and 


*» 


dv  dja(v) 


(2) 


is  the  well-known  dielectric  function  for  longitudinal  plasma  oscillations,  with 
7*(v)  representing  the  equilibrium  velocity  distribution  function  of  species  <r. 


Here  and  henceforth,  the  bar  and  the  tilde  indicate  equilibrium  quantities  and 
perturbations,  respectively.  For  a  given  value  of  q,  Eqs.  (1)  and  (2)  as  they  stand 
are  valid  only  in  the  upper  u  half-plane;  if  Imu  <0,  the  path  of  integration 
must  usually  be  distorted  in  the  complex  o-plane,  so  as  to  ensure  proper  analytic 
continuation  with  respect  to  u. 

Inverting  V{q,u)  into  V(q,t)  one  finds  that  the  latter  contains,  among 
other  possible  contributions,  terms  proportional  to  exp(— tw„t),  where  the  mode 
frequencies  u„  (x/  1,2, ...)  are  the  zeros  of  the  dielectric  function, 

D(q,u„)  —  0,  (3) 

and  hence  also  represent  poles  of  V’fa.w)  as  given  by  (1). 

Ever  since  its  appearance  in  the  literature,  the  infinite-plasma  dispersion 
relation  (3)  (or  some  analogous  dispersion  relation  for  more  complex  veloc¬ 
ity  structures)  has  been  used  by  numerous  authors  for  investigating  the  wave- 
propagation  and  stability  properties  of  a  great  variety  of  plasma  configura¬ 
tions.*'*  As  long  as  the  wavelengths  in  question  are  much  shorter  than  the 
relevant  macroscopic  dimensions,  this  approach  is  usually  believed  to  be  jus¬ 
tified.  However,  due  to  its  relative  simplicity,  and  owing  to  the  fact  that  a 
proper  treatment  of  bounded  plasma  systems  is  usually  much  more  difficult  and, 
hence,  not  readily  feasible  (the  eigenmodes  are  no  longer  identical  with  single 
Fourier  modes),  Landau's  infinite-plasma  result  is  often  also  applied  to  situa¬ 
tions  where  a  bounded-system  analysis  would  be  more  appropriate.4*'  (Rognlien 
et  al.r  found  that  the  eigenmodes  of  systems  described  by  fiuid-type  equations 
can  under  certain  circumstances,  but  not  generally,  be  constructed  from  the 
roots  of  the  corresponding  infinite-system  dispersion  relation.  At  the  same 
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time,  however,  they  pointed  out  that  their  treatment  did  not  in  general  apply 
to  collisionless  plasmas  described  by  kinetic  equations.) 

Thus,  the  question  we  wish  to  answer  in  the  present  work  is  the  following; 
What  is  the  equivalent  of  Landau's  fundamental  result  (1)  for  linear  pertur¬ 
bations  in  collisionless  bounded  plasma  systems,  where  not  only  equilibrium 
plasma  properties  and  initial  conditions,  but  also  boundaries  and  external  cir¬ 
cuits  have  to  be  taken  into  account?  May  we  expect  explicit  results  similar  to 
(1)  (which  would  be  the  most  favorable  case),  or  will  we  end  up  with  some¬ 
thing  considerably  more  complex?  If  the  latter  is  true,  how  can  we  still  extract 
information  on  dispersion  and  stability  properties? 

We  decide  to  solve  this  problem  for  a  system  that  is  sufficiently  simple  to 
allow  for  a  basically  analytic  treatment,  but  still  complex  enough  to  account 
for  typical  bounded-system  effects. 

The  system  we  choose  is  the  collisionless,  one-dimensional  plasma  diode 
with  “thin"  sheaths  (in  a  sense  to  be  specified  in  Subsec.  II.A)  but  otherwise 
allowing  for  very  general  equilibrium,  boundary,  and  external-circuit  condi¬ 
tions.  The  major  geometrical  assumptions,  namely  one-dimensionality  and  thin 
sheaths,  are  not  believed  to  be  crucial  to  the  whole  (integral-equation)  approach 
as  such  (cf.  Subsec.  Dd.A)  but  are  certainly  helpful  in  this  first  stage  of  devel¬ 
opment  of  the  theory. 

The  present  work  is  divided  into  two  parts.  In  Part  I  (this  paper),  the 
general  formalism  is  developed  (Sec.  II)  and  conclusions  of  a  general  nature  are 
drawn  (Secs.  01,  IV).  Part  II*  presents  a  detailed  and  quantitative  application 
of  the  formalism  to  a  very  fundamental  special  case,  namely  to  an  extended 
Pierce-type  problem  involving  a  noo-trivial  external  circuit.  This  is  a  striking 
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example  of  how  a  plasma  that  is  stable  in  the  infinite-system  limit  may  become 
violently  unstable  if  part  of  a  bounded  system.  The  numerical  results  presented 
in  Part  II  also  demonstrate  that  the  external  circuit  may  have  drastic  effects 
on  the  stability  and  dispersion  properties  of  the  whole  system,  a  fact  which  is 
frequently  ignored. 

The  rest  of  this  paper  (Part  I)  is  organized  as  follows.  Subsection  II.A  de¬ 
scribes  the  model  considered  and  introduces  the  basic  equations  (Vlasov-Poisson) 
for  the  plasma  region.  In  Subsec.  II.B  these  equations  are  Laplace  transformed 
with  respect  to  t  and  x,  the  perturbations  of  the  velocity  distribution  functions 
are  eliminated,  and  the  resulting  form  of  Poisson’s  equation  is  Laplace  inverted 
with  respect  to  2.  In  Subsection  II.C,  the  condition  of  total-current  conservation 
is  used  to  “link”  the  plasma  to  the  external  circuit  at  the  left-hand  electrode, 
and  Poisson's  equation  is  cast  into  a  form  symmetric  with  respect  to  both  elec¬ 
trodes.  Subsection  ILD  presents  a  very  general  class  of  boundary  conditions, 
whence  a  set  of  integral  equations  is  derived  for  the  perturbations  of  the  veloc¬ 
ity  distribution  functions  at  the  plasma-sheath  boundaries.  In  Subsec.  II.E  we 
adopt  a  very  general  class  of  external-circuit  conditions,  thus  closing  the  system 
of  equations  needed  to  determine  the  dynamics  of  our  linearly  perturbed  diode 
system. 

Subsection  Q.F  reviews  the  structure  of  the  problem,  which  in  the  most 
general  case  involves  2  +  2n„  coupled  integral  equations  in  z  and  v  for  the  fol¬ 
lowing  quantities,  which  are  the  time  Laplace  transforms  of  the  corresponding 
physical  perturbations:  y,(w)  (external-circuit  current  density),  E(x,  w)  (electro¬ 
static  field),  /f(»  >  0,w)  and  ~/*{v  <  O.w)  (velocity  distribution  functions  for 
the  plasma-bound  particles  at  the  left-hand  and  right-hand  plasma  boundaries. 
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respectively).  These  integral  equations  basically  represent  the  bounded-system 
analogue  to  Landau’s  infinite-plasma  result  (1)  and  can  be  easily  shown  to  re¬ 
duce  to  the  latter  in  the  appropriate  limit.  Explicit  solutions  for  the  time 
Laplace  transforms  jt(u )  etc.  may  be  expected  only  for  '‘sufficiently’’  simple 
systems.  However,  the  intrinsic  eigenfrequencies  w„  (as  opposed  to  frequencies 
generated  externally)  are  shown  to  be  the  zeros  of  the  coefficient  determinant  of 
the  up  to  2+2na  integral  equations  in  a  suitable  basis-set  representation.  Thus, 
there  is  still  a  systematic  way  of  analyzing  stability  and  dispersion  properties 
for  any  given  application. 

The  discussion  in  Subsec.  IH.A  leads  to  the  main  conclusion  that  the  integ¬ 
ral-equation  approach  adopted  here  is  likely  to  be  extendible  to  more  general 
and  complex  problems.  In  Subsec.  m.B  we  review  some  pertinent  literature  and 
hence  conclude  that  the  present  method  includes  a  significant  class  of  previous 
treatments  as  special  cases  and,  in  addition,  has  potential  for  investigating 
problems  that  have  hardly  been  touched  in  the  previous  literature.  Finally, 
Sec.  IV  summarizes  our  main  conclusions. 


....  Sft-gfc 


IL  METHOD 


A.  Model  and  basic  equations 

We  consider  a  one-dimensional  diode  as  shown  in  Fig.  1.  The  surfaces  of 
the  (ideally  conducting)  electrodes  are  located  at  z  0  (“left-hand  electrode”) 
and  z  ™  L  (“right-hand  electrode").  Let  the  intervening  space  be  filled  with  a 
collisionless  plasma  consisting  of  n„  particle  species,  and  let  the  far  ends  of  the 
electrodes  be  connected  through  an  external  circuit  with  specified  properties 
(Subsec.  n.E). 

For  the  d.c.  state  to  be  perturbed  (henceforth  referred  to  as  “equilibrium" ) 
we  assume  a  uniform  plasma  with  a  constant  plasma  potential  V r,  extending 
over  the  region  zi  <  z  <  zr.  In  the  regions  0  <  z  <  zi  and  z,  <  z  <  L 
we  allow  for  space-charge  sheaths,  which  are  crucial  in  shaping  the  equilibrium 
velocity  distribution  functions.*'1*  With  respect  to  the  perturbations,  these 
sheaths  are  assumed  to  be  “thin"  (*<  —  0,z,  -►  L).1*  This  approximation  is 
justified  if  (i)  the  typical  times  of  particle  transit  through  the  sheath  regions 
are  short  compared  with  the  characteristic  time  scales  of  the  perturbations  (cf. 
Subsec.  HD),  and  (ii)  the  sheath  widths  are  small  with  respect  to  the  typical 
scale  lengths  of  the  perturbations  (cf.  Subsec.  HE).  Globally  speaking,  these 
requirements  restrict  the  validity  of  our  analysis  to  perturbations  which  do  not 
crucially  depend  on  the  fine  structures  of  the  sheaths,  and  for  which  the  sheath 
impedances  can  be  neglected. 

The  thin-sheath  approximation  is  of  relevance,  e.g.,  for  longitudinal  modes 
in  a  single-ended  Q  machine14,14  at  “moderate"  interelectrode  biases,  whereas 
for  “very  high"  biases  the  sheath  widths  may  no  longer  be  negligible.1*17  As  an 
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example,  Fig.  1  shows  a  “one-emitter"  potential  distribution,  which  represents 
one  of  the  self-consistent  equilibrium  configurations  of  the  ideal  one-emitter 
plasma  diode  or  single-ended  Q  machine.*'11 


In  the  plasma  region,  the  small-amplitude  longitudinal  oscillations  we  wish 
to  study  are  described  by  the  linearized  Vlasov  and  Poisson  equations: 


*21 

at 

dE 

—  s 

Qz 


+.£—4  r.E 

dx  m°  * 
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4*  51  *  I  dv  f  ' 
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where  E(x,t)  =*  E(z,t)  is  the  electrostatic  field,  ma  and  Ja(z,  v,t)  =  /^(v)  + 

•  ff 

f  ( z ,  v,  f)  are  the  particle  mass  and  velocity  distribution  function  of  species  <r, 
and  d]a(o}/do. 


Equation  (4)  has  the  following  formal  solution,  which  will  be  needed  later 
on: 

[  <**£(*•*)  (#) 

m*  v 

where  {*  *  x  +  tr(t  —  t),  v  *  »,t},  with  t  as  the  basic  parameter,  is  any 
point  along  the  equilibrium  particle  trajectory  passing  through  the  given  point 
(a,  v,  t),  and  {x0  =  z  +  v(t0  —  t),  v0  =  »,  t0 }  is  the  starting  point  for  integration 
along  the  trajectory. 

The  problem  to  be  solved  consists  in  predicting  the  linear  evolution  in  time 
of  the  perturbed  diode  system  for  given  equilibrium,  initial,  boundary,  and 
external-circuit  conditions.  In  the  rest  of  this  section  we  present  a  very  general 
method  for  tackling  this  problem,  which  implies  introducing  several  additional 
relations  and  unknowns.  The  following  subsections  describe  the  single  steps 


-10- 


which  will  eventually  (Subsec.  ILF)  lead  to  a  set  of  coupled  integral  equations 
for  the  time  Laplace  transforms  of  the  dynamic  variables  chosen. 

B.  Laplace  transformations  and  x-inversion 

In  order  to  solve  the  linearized  basic  equations  in  the  plasma  region,  we 
choose  to  make  use  of  Laplace  transformations  in  both  space  and  time.  Let 
these  be  defined  by 

5(f)-  /  dze^S(z),  5(a)-  /  S(?)  (7a,b) 

J 0  J— oo+t(!inf)i  ** 

and 

rfoo+Wmvh  fa  . 

<HeMT{t),  T(«)-  /  •—  e“,w<  T(w) ,  (8a,b) 

J—c o+((Imw)i  ** 

where,  as  usual,  Im?,  (Im?)i,  Imu,  and  (Imw)i  must  be  “sufficiently"  positive. 
A  major  advantage  of  using  these  transformations  lies  it.  the  fact  that  they 
permit  one  to  handle  boundary  and  initial  conditions  in  a  very  natural  and 
convenient  manner.  Note  that,  although  the  spatial  transformation  is  formally 
applied  to  the  whole  half-space  (0  <  z  <  oo),  only  the  plasma  region  (zt  < 
z  $  *r)  will  ultimately  be  of  physical  interest. 

To  our  knowledge,  the  only  previous  solution  of  the  linearized  Vlasov  equa¬ 
tion  by  means  of  the  double  Laplace  transform  method  has  been  given  by 
Evans.1*  In  Refs.  19  and  20  this  method  was  also  applied  in  the  context  of 
cold-fiuid  type  systems.  However,  all  of  these  treatments  were  concerned  with 
semi-infinite  systems,  which  may  be  viewed  as  special  limiting  cases  of  the  diode 
configuration  considered  here.  A  preliminary  version  of  the  present  method  was 
given  in  Ref.  21. 


so  that  (11)  can  b«  easily  re-written  in  terms  of  V[q,u).  Let  us,  just  for  the 
present  purpose  and  without  lack  of  generality,  transform  the  interval  0  <  x  <  L 
into  the  new  one  — L/2  <  *  <  Lj 2,  with  L  —  oo.  Since  Landau1  used  Fourier 
integrals  to  represent  his  perturbations,  the  latter  were  tacitly  assumed  to  be 

-  (J 

localized  in  space.  In  particular,  this  is  the  case  for  /  (z,  v,  t),  E(z,  t ),  and 
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V{z,t),  so  that  these  quantities  all  become  zero  at  x  ■  ±L/2  as  L  tends  to 
infinity.  This  means  that  the  quantities  £;(<•>)  and  7*(v,u)  in  Eq.  (11)  must 
vanish.  Finally,  letting  Imq  become  real  leads  us,  via  (13),  to  Eq.  (1). 

By  inspection  of  Eq.  (11)  we  see  that  the  dependence  of  E{f,u)  on  q  is 
explicit  once  the  equilibrium  distribution  functions  ?*{v)  and  the  initial  per* 
turbations  ft(z,  *)  have  been  specified.  Hence,  we  can  now  apply  to  (11)  the 
inverse  spatial  transformation  (7b),  which  leads  to 

«<x>  J 

E(x,«)-«~4#53  s' /  — *'(*,*, w)  +  »*j(*,w)E((w) 

..  (14) 

“4  r£e"/  dvkt(z,a,u)f*(v,u). 

ml 

The  known  functions  it*,  kj  and  kt  introduced  here  are  defined  in  Eqs.  (Al) 
through  (A3)  of  Appendix  A.  There  we  have  also  listed  a  number  of  other 
functions  which  will  be  encountered  in  the  course  of  this  study  and  may  be 
considered  as  being  known,  in  the  sense  that  they  can,  in  principle,  be  evaluated 
once  the  equilibrium,  initial,  boundary,  and  external-circuit  conditions  have 
been  specified. 

In  deriving  Eq.  (14)  we  have  interchanged  the  order  of  performing  the 
inverse  spatial  Laplace  transformation  (7b)  and  the  velocity  integration.  It  is 
not  difficult  to  justify  this  procedure  for  the  case  of  a  finite  number  of  cold 
beams.  Moreover,  since  any  plasma  may  be  approximated  by  a  sufficiently 
large  number  of  cold  beams,  we  conclude  that  it  is  also  admissible  for  plasmas 
with  continuous  velocity  distribution  functions. 

If  the  perturbations  2?<(w)  and  7*[v,u)  were  explicitly  known  as  functions 
of  v,  Eq.  (14)  could  now  be  Laplace  inverted  for  E(z,t),  so  that  our  problem 
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would  be  paractically  solved.  This  is,  ia  fact,  possible  at  least  for  Landau's 
infinite-plasma  problem  (cf.  Sec.  I),  but  in  general  the  above  functions  must 
be  eliminated  or  determined  via  suitable  additional  relations.  In  the  following 
subsections,  these  relations  are  introduced  and  incorporated  into  our  formalism. 

C.  Total-current  conservation  and  symmetric  form  of  Poisson’s 
equation 


In  a  one-dimensional  system  of  the  form  considered  here,  the  total  (i.e., 
convection  plus  displacement)  current  density  depends  on  time  only.22  At  the 
left-hand  electrode  this  condition  assumes  the  form 


(15) 


4r  dt 

where  jt(t)  is  the  perturbation  of  the  external-current  density,  and  /,(f)  is  the 
perturbation  of  the  convection  or  particle  current  density  at  x  =**  zt: 


J‘l(*)  *  51  /  dvvf"(v,t). 

J—OO 


(18) 


Inserting  (16)  into  (15)  and  applying  the  time  Laplace  transformation  (8a) 

we  find 


■  _  »»  roo 

Ei(u)  —  -  Eh  +  4wjt(u)  -  4r  t°  /  ivvj^v,  w] 


(17) 


The  constant  Eh  s  E(z  —  zi,  t  —  0)  can  be  calculated  in  terms  of  other  initial 
values  by  integrating  (5)  twice  with  respect  to  z,  at  t  =  0: 


E,i~-~gi-j-  £  Ax  j‘  dz'  r  */?(*'.») 


(18) 


where  AV^*  5  AK,(t  —  0),  A V,(t)  being  the  perturbation  of  the  potential 
drop  across  the  plasma  region  (cf.  Fig.  1  and  Eq.  (12)): 


> - _ 
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izE(z,t). 


According  to  Sec.  II.E  below,  the  initial  potential-difference  perturbation  AV f{ 
is  essentially  determined  by  the  initial  state  of  the  external  circuit. 

We  now  insert  Ei(u)  as  given  by  Eq.  (17)jnto  Poisson’s  equation  (14)  and, 
after  rearranging  terms,  obtain 

E(x,u)  »*4(z,w)  +  *g(z,w)j,(w) 


*»  -OO  .  g 

+  4r  ^  t9  I  doki(z,v,u)fl(v,u). 


where  the  (known)  functions  kt,  kg,  and  k«  are  defined  by  Eqs.  (A4)  through 
(A6).  Although  Eq.(20)  contains  the  same  number  of  unknowns  as  Eq.(14), 
progress  has  been  made  in  that  Ei(u)  has  been  replaced  with  the  more  useful 
quantity  which  establishes  a  link  between  the  plasma  and  the  external 

circuit. 

Equation  (20)  contains  a  velocity  integral  to  be  evaluated  sti»  zj,  but 
no  such  integral  for  z  *«  zr.  It  is  thus  asymmetric  with  respect  to  the  two 
electrodes,  none  of  which  is,  however,  a  priori  assigned  a  preferred  role.  We 
choose  to  remove  this  asymmetry  by  expressing  (v  <  0,  w)  in  terms  of  f°(v  < 
0,u),  to  which  end  we  first  evaluate  (6)  at  z  ®  zj  —  0: 

/f(»  CO,*)-  U(t  +  £)  fdiE{z,t  +  f )} 

»  l  v  m  v  JL  v  ) 

+uH  -  i)  {/'(-«. .)  -  iSjd  +  i>}  . 

where  U  is  the  Heaviside  unit  step  function.  Obviously,  the  term  with  U(t+L/v) 
represents  the  contribution  of  those  particles  which  took  off  at  the  right-hand 
electrode  at  a  time  t«  >  0.  The  second  term  describes  the  influence  of  those 


(slower)  particles  which  at  time  i9  »  0  were  already  on  their  way  towards 
the  left-hand  electrode  and  thus  are  the  immediate  carriers  of  the  information 
contained  in  the  initial  conditions.  The  time  Laplace  transform  (8a)  of  Eq.  (21) 


/,(»  <  o,w)  —  -  i*7(t>,w)+ 

v  rn  v 

+  exp(-»w-)  /'(e.w). 

V 


exp(—iu—)E(x',u) 


(22) 


For  later  use  we  also  write  down  the  analogous  relations  for  the  right-hand 
electrode,  namely 


r 


U*  >0,0  —  -  t )  |/?(L  £  *  HEd.i)  J 


(23) 


and  its  time  Laplace  transform 


f’r(v  >  0,w)» iexp(»«-)*^iu,w) 

0  V 

—  exp(iw-)  f  dx'  exp{—iu  —  )E[x',u)  (24) 

m  v  v  J  0  v 

+  exp(w^)/f  (»,«). 

V 


Inserting  (22)  into  (20)  we  obtain 

E(x,u)  *  kt(x,w)  +  ki{x,v)j,{u)  +  J*  dx'k,(x,x\v)E(x',u) 

+  4*  J~]  *{£  dvkt{X'V’  v)/Uv<u)  (25) 

+  J  dvk,{x,v,u)  exp(-iw£)/"(«;,w)j  , 


with  the  newly  introduced  known  functions  k9  and  kt  defined  in  Appendix  A. 
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Equation  (25)  is  the  more  symmetric  version  of  Poisson's  equation  that 
we  have  been  looking  for.  It  contains  two  source  terms  corresponding  to  the 
particles  entering  the  plasma  region  at  the  plasma-sheath  boundaries.  If  the 
boundary  perturbations  /f(o  >  0,  w)  and  /*(v  <  0,  w)  were  explicitly  known, 
(25)  would  represent  a  Fredholm-type  integral  equation  for  E(x,u)  only,  with 
w  and  j'4(w)  as  parameters.  However,  in  many  cases  of  interest  these  functions 
will  not  be  known  a  priori  but  have  to  be  determined  with  the  help  of  the 
relevant  boundary  conditions.  Hence,  the  latter  must  now  be  specified.  In  the 
following  subsection  we  introduce  a  very  general  class  of  boundary  conditions 
and  therewith  derive  additional  equations  for  our  unknowns. 

D.  Boundary  conditions  on  particles  and  resulting  integral  equa¬ 
tions 


Let  us  assume  that,  at  the  plasma-sheath  boundaries  z  »  zi  and  z  ■  zr, 
the  perturbations  of  the  velocity  distribution  functions  of  the  plasma^bound 
particles  are  related  to  those  of  the  sheath-bound  particles  as  follows: 


7f(»  >  0,t)»  /^M)+  £  /  do' bf*'(v,v') /f  (v\ t) 

<7*— 1  J-OO 


(26a) 


fUv  <  0,*)  — 7^M)+  X)  [  dv'  b“'(v,v')  /"(»’,*), 


(28b) 


o'— 1 


where  f*t  and  are  externally  generated  (and,  hence,  explicitly  given)  per¬ 
turbations,  and  the  functions  b**  (»  >  0,  v'  <  0)  and  6**  (»  <  0,  v1  >  0) 
essentially  represent  the  probabilities  for  a  sheath-bound  particle  of  species  o’ 
with  velocity  o'  to  “produce"  a  plasma-bound  particle  of  species  a  with  velocity 
o.  The  time  Laplace  transforms  of  (26a,b)  are 


A 


As  an  example  of  an  external  perturbation  mechanism  let  us  mention  the 
fluctuations  in  the  flux  of  neutral  alkali  atoms  incident  on  a  hot  metal  electrode 
in  a  Q  machine,  which  may  cause  variations  in  the  flux  of  the  plasma-bound 
ions  produced  by  contact  ionization. 

The  homogeneous  terms  in  (26a,b)  are  capable  of  describing  a  wide  class 
of  generalized  reflection  processes  that  a  sheath-bound  particle  of  species  o’ 
incident  with  velocity  o'  can  undergo.  Among  them  are  some  simple  ones  most 
commonly  used  (cf.  Subsec.  m.B).  If,  e.g.,  the  sheath-bound  particle  is  always 


absorbed  without  releasing  any  other  particle,  we  have 

bl'\v  >  0,  o'  <  0)  —  0  (28a) 

b"’(o  <  0,o' >  0)  — 0  (28b) 

for  all  permitted  values  of  9  and  o.  For  specular  reflection,  the  b's  become 

6f*'(o  >  0,  o'  <  0)  -  6av.  S(v  +  o')  (29a) 

b’/(v  <  0,  o'  >  0)  »  6ov.  6(o  +  o') ,  (29b) 


where  S/y  is  the  Kronecker  delta  symbol  and  6(x)  is  the  Dirac  delta  function. 
Clearly,  the  boundary  conditions  (26)  can  also  describe  processes  like  inelastic 
reflection  or  secondary-particle  emission,  but  these  should  be  discussed  in  spe¬ 
cific  future  applications  rather  than  in  the  present  general  outline  of  the  theory. 
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la  (28),  the  perturbations  of  the  velocity  distribution  functions  for  the 
plasma- bound  particles  are  evaluated  at  the  same  time  t  as  those  for  the  sheath- 
bound  particles.  This  involves  two  major  physical  assumptions,  namely  that 
both  (i)  the  time  needed  for  a  particle  to  traverse  the  sheath  regions,  and  (ii)  the 
time  scales  of  the  generalized  reflection  processes  taking  place  at  the  electrode 
surfaces  proper,  are  negligibly  small  in  comparison  with  the  time  scales  of  the 
dynamic  phenomena  of  interest.  However,  these  ‘instant-reflection’'  conditions 
are  likely  to  be  satisfied  in  many  cases  of  practical  interest. 

It  should  be  borne  in  mind  that  the  "boundary”  conditions  (28)  are  not 
only  determined  by  the  electrode  surface  properties  but  also  account,  in  a  global 
manner,  for  the  sheath  regions.  (This  global  treatment  is  made  possible  by 
the  assumption  in  Subsec.  II.A  of  "thin”  sheaths,  whereas  a  finite-sheath  the¬ 
ory  would  require  a  more  refined  approximation  to  the  solution  of  the  Vlasov- 
Poisson  system  in  the  sheath  regions.)  Consider,  e.g.,  a  singly  charged  positive 
ion  of  mass  m  moving  to  the  right  in  the  configuration  shown  in  Fig.  1.  If 
o  <  (2 e(V l  —  V ,) /m\1/*,  the  ion  is  reflected  somewhere  in  the  right-hand 
sheath  region  (xi  <  z  <  L),  which  in  terms  of  (28)  counts  as  a  specular  reflec¬ 
tion,  cf.  (29b).  In  the  opposite  case  the  ion  hits  the  electrode  surface  proper  and 
there  may  be  neutralized,  undergo  inelastic  reflection,  or  even  release  secondary 
particles.  Clearly,  the  description  of  these  latter  processes  will  require  a  form 
of  b*a'  considerably  more  complex  than  (29b). 

We  are  now  in  a  position  to  derive  the  additional  relations  required  for 

m  0  •  y 

determining  the  functions  /,  (o  >  0,w)  and  fr(v  <  0,«),  which  occur  in 
Poisson's  equation  (28).  Let  us  start  with  fj(v  >  0,w)  as  given  by  (27a),  and 

m  0  -  y 

express  fx(v  <  0,w)  on  the  r.h.s.  step  by  step  in  terms  of  Jf{v  <  0tw), 


*  -V. 


( 
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f,(v  >  0, w),  and  fi(v  >  0,  w)  by  successive  application  of  Eqs.  (22),  (27b), 
and  (24).  What  comes  out  is  the  integral  equation 

rt 


f"(v  >  0,w)  *  k"9(v,u)  +  J  dx'  kf^x'.u.w)  £(*',«) 


(T'—l 


(30) 


)• 


In  an  analogous  manner  we  can  express  fr(v  <  0,u)  consecutively  in  terms 
of  /"{v  >  0,  w),  f*(*  >  0,w),  }1(v  <  0,w),  and  ~f*(v  <  0,  w)  by  using  Eqs. 
(27b),  (24),  (27a),  and  (22).  The  result  is 


far(v  <  0,w)  —  *",(»,«)  +  J  dx'  *?4(x',»,w)  E(z',u) 


+  £  f  dv' k*Z\v,v',U)  /;'(»',  w). 


(31) 


Again,  the  newly  introduced  known  functions  are  defined  in  Appendix  A. 


Equations  (25),  (30),  and  (31)  constitute  a  system  of  l+2n,  coupled  integral 
equations  in  x  and  v  for  the  1  +  2 n„  functions  E(z,u),  f\(v  >  0,w),  and 
/,(»  <  0,  w).  Since  both  w  and  j,(u)  enter  as  parameters,  yet  another  relation 
is  required  in  order  to  specify  the  w-dependence  of  the  above  functions  and  to 
thus  make  them  invertible  into  E(x,t)  etc.  In  the  following  subsection  we  show 
how  this  relation  can  be  obtained  by  combining  the  overall  potential  balance 
with  the  external-circuit  condition. 

E.  Overall  potential  balance  and  external-circuit  condition 

So  far  we  have  essentially  concentrated  on  the  plasma  region,  with  y4(w) 
entering  as  a  free  parameter,  cf.  Eq.  (25).  Thus,  in  order  to  close  the  prob¬ 
lem  we  need  yet  another  suitable  equation  that  will  completely  specify  the  in¬ 
dependence  of  the  quantities  involved.  In  the  present  subsection  we  derive  this 
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equation  by  combining  (i)  the  overall  potential  balance  in  the  diode/external- 
circuit  system,  (ii)  the  potential  balance  in  the  external  circuit,  and  (iii)  the 
condition  of  zero  sheath  impedance,  which  is  consistent  with  the  thin-sheath 
approximation. 

The  overall  balance  for  the  potential  perturbations  may  be  written 

AK.(t)  -  AV-#1(<)  +  A7,(«)  +  A7„(t) ,  (32) 

where  A7„(*)  -  7(*,,*)  -  7(0,*)  and  AV„(*)  -  V(L,t)  -  V(zr,t)  are  the 
perturbations  of  the  potential  drops  across  the  left-hand  and  the  right-hand 
sheath  regions,  respectively.  The  potential  balance  for  the  external  circuit 
(“external-circuit  condition")  we  choose  to  write  in  the  form 

;*,(»),  (33) 

where  Z(t)  is  a  general  linear  impedance  operator  that  may  be  considered  known 
once  the  external-circuit  properties  have  been  specified. 

In  order  for  (32)  to  be  useful  in  the  context  of  our  preceding  considerations 
(where  the  basic  equations  were  solved  for  the  uniform  plasma  region),  AV,i(*) 
and  ■3V’tf(t)  must  be  calculated.  In  a  finite-sheath  theory  this  would  imply  at 
least  approximate  integration  of  the  Vlasov-Poisson  system  in  the  (non-uniform) 
sheath  regions.  However,  in  the  thin-sheath  approximation  considered  through¬ 
out  this  paper  (cf.  Subsec.  II.A)  we  may,  by  definition,  set 

A7.,(t)-0-  AV-„(t),  (34) 

so  that  (32)  reduces  to 

A7.(*)-  A7,(»).  (35) 

This  means  that  the  perturbations  of  the  potential  drops  developing  over  the 
thin  sheath  regions  are  neglected  against  those  building  up  over  the  much  more 


extended  plasma  region.  (This  behavior  is  exactly  opposed  to  equilibrium  con¬ 
ditions,  where  the  main  potential  drops  occur  over  the  sheath  regions.) 

The  operator  Z  introduced  in  Eq.  (32)  is  assumed  to  be  linear  but  otherwise 
completely  general.  However,  for  many  (if  not  all)  cases  of  practical  interest  Z 
will  be  of  such  a  form  that  the  time  Laplace  transform  of  (33)  becomes 

AV.(«)  =  i^,(w)  +  Z>)J>»).  (36) 

This  class  of  external  circuits  (which  is  still  very  general  in  that  it  includes 
the  usual  RLC  circuits  with  a.c.  sources)  is  the  one  considered  henceforth.  It 
may  be  useful  to  note  that  /^(w)  reflects  the  effect  of  initial  conditions  and/or 
external  sources,  whereas  Zt(u)  is  a  generalized  impedance.  The  initial  state  of 
the  external  circuit  also  specifies  the  initial  potential-difference  perturbations 
&Vti  and  AVfi. 

An  explicit  treatment  of  an  exemplary  external  circuit  is  given  in  Part  II,* 
where  the  above  aspects  are  worked  out  in  detail. 

Equation  (36),  with  AV,(w)  *»  AV,(w),  closes  the  system  of  equations 
governing  our  perturbational  problem,  whose  general  structure  is  discussed  in 
the  following  subsection. 


F.  General  structure  of  the  problem  and  basis-set  expansion 

Equations  (30)  (combined  with  (35)  and  (IS)),  (25),  (30),  and  (31)  form  a 


closed  set  of  2+2n»  coupled  integral  equations,  from  which  the  2+2n*  unknown 
time  Laplace  transforms  jt(v),  E(*,u),  /t  (o  >  0,  w),  and  fr(v  <  0 ,u)  can, 
in  principle,  be  determined  as  functions  of  x,  v,  and  u.  (Once  these  solutions 
have  been  established,  /(x,  o,u)  could  be  calculated  via  (S).)  The  last  step 
in  solving  the  complete  perturbational  problem  then  consists  in  applying  the 
inverse  time  Laplace  transform  (8b)  to  the  above  functions,  which  will  yield 
the  time-dependent  perturbations  y4(t),  E(z,t),  f*(v  >  0,t),  f"(v  <  0,t),  and 
finally,  via  (6),  /'(x,  v,  t)  (<r  —  1, . . . ,  n,). 

In  order  to  reduce  the  problem  to  its  basic  structure,  we  re-write  the  above 
equations  more  concisely  as  follows: 

Z.{»)  j»  +  jf  *  «fa'  E(z\  w)  -  -V-*(u/)  (37) 

-  *i(*, «)  i.(«)  +  E(x ,  u)  +  S®(x,  (x'J,  u)  E([x'],  u)  ’ 

+  5Z  l»  >  °!> ")  n (W. ») 

+  5Z  <  °]>")  —  *•(*»«) 


Sf(»  >  0,(x'J,«)  E{[z'],u) 


+/f(*,u)+  ^2  vr'(u,K>0|,W)/r(Kj,«)-*r.(.,«) 


(30) 


S*(»  <  0,  (x'|,w)  E([x'],u) 

+  /'(*,  u)  +  £  l«'  <  0],  w)  /'  ([•'!,«)  W) 


(40) 


ittiiii  ii  - 1  n  ii  i  ^ni 


-M- 


where  S«,  S j,  S,  are  known  space-integral  operators  and  "Vjj,  "VJ,.,  \f<r\ 
are  known  velocity-integral  operators.  These  operators  are  defined  in  Eqs.  (A16) 
through  (A22)  of  Appendix  A.  The  square  brackets  indicate  the  variables  of 
integration. 

The  integral  equations  (37)— (-40)  may  be  viewed  as  the  bounded-system 
analogue  to  the  infinite-system  result  (1),  which  answers  the  question  raised 
in  Sec.I.  It  thus  turns  out  that  in  general  (i.e.,  for  sufficiently  non-degenerate 
cases)  it  is  no  longer  possible  to  obtain  the  time  Laplace  transforms  of  the 
perturbations  explicitly  as  functions  of  u.  However,  in  Subsec.  II.G  it  will  turn 
out  that  there  still  exists  a  systematic  way  of  analyzing  stability  and  dispersion 
properties. 

Depending  on  the  degeneracy  of  a  specific  configuration  considered,  some  of 
the  equations  (37)-{40)  may  be  trivially  satisfied  (as  is  the  case  in  the  extended 
Pierce-type  problem  treated  in  Part  II8),  so  that  only  an  appropriate  subset 
thereof  must  be  actually  considered. 

In  practical  applications,  when  the  coefficient  functions  and  integral  oper¬ 
ators  occurring  in  Eqs.  (37)-(40)  must  be  explicitly  established,  it  is  convenient 
to  distinguish  between  four  levels  of  quantities  as  illustrated  by  Table  I.  The 
level-one  quantities,  which  determine  the  equilibrium,  initial,  boundary,  and 
external-circuit  conditions,  must  be  specified,  while  any  higher-level  quantity 
can  be  constructed  from  lower-level  ones.  For  each  quantity  is  indicated  in 
Table  I  an  equation  or  section  where  it  is  defined  or  introduced.  The  level-four 
quantities  are  the  functions  and  operators  occurring  in  Eqs.  (37)-(40). 

Once  established  explicitly,  Eqs.  (37)-(40)  will  in  many  cases  be  too  com¬ 
plicated  to  allow  for  closed  analytical  solutions.  It  may  then  be  convenient  to 
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expand  the  functions  and  operators  involved  in  suitable  sets  of  basis  functions. 
This  permits  one  to  transform  the  integral  equations  into  a  system  of  linear  al¬ 
gebraic  equations  for  the  expansion  coefficients  of  the  unknown  functions.  The 
coordinate  ranges  to  be  covered  are  0  <  v  <  o-,,,  and  v*,,*  <  v  <  0,  where 
the  cutoff  velocities  and  om<)(  must  be  chosen  so  as  to  cover  all  phenom¬ 
ena  of  interest  in  velocity  space.  Let  there  be  given  three  denumerably  infinite, 
complete  sets  of  linearly  independent  (but  not  necessarily  orthogonal  and/or 
normalized)  basis  functions  {pK(x)),  >  0)},  and  {p~[v  <  0)},  defined 


in  the  above  intervals,  respectively,  with  *,  X,  (i  ■*  ...1,2,... 

involved  may  then  be  expanded  in  one  of  the  forms 

.  Any  function 

f(z)  •>  /(*)  v»«(i) 

(41a) 

/,(«  >  0)  -  /,(X) 

(41b) 

fT(0  <  0)  —  <p~(v) 

(41c) 

or  straightforward  extensions  thereof,  e.g. 

g{x,v  <  0)  —  pjx)  p~{v). 

(42) 

Note  that  we  have  adopted  the  Einstein  convention  (double  indeces  imply  sum¬ 
mation)  for  *,  X,  and  p  (but  not  for  the  species  index  v).  .\s  usual,  these 
expansions  permit  us  to  represent  continuous  functions  and  operators  in  terms 
of  their  respective  coefficient  vectors  or  matrices.  Generally,  the  latter  are  in¬ 
finite  but  in  practice  must  always  be  truncated,  the  number  of  terms  retained 
usually  being  a  compromise  between  accuracy  requirements  and  technical  lim¬ 
itations. 


A 
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where  0,  <&,  kg,  kg,  k*0,  kfj,  E,  ff,  t°  are  infinite  column  (subsectors  with 
components  0,  L~l  dxifijz),  kg(*,w),  Ag(/e,w),  fe[0(\,v),  *u(m.w).  E{k,u), 
/"(X.w),  /'(/»,  w),  respectively,  the  superscript  T  indicates  transposed  vectors, 
and  0,  L  £o,  §j<  S.t<  YZi,  Zo,.  Kfr'<  Yf/'  are  inflate  (sub)matrices  with 
elements 0.  £««/,  S0(k,k',u),  Sf(X,ie',w) «',«),  V^(k, X',w),  V0*(k,/*',w), 
V"  (X,  X'.u),  respectively.  Let  us  recall  that  in  practice  these 

vectors  and  matrices  must  always  be  truncated. 

As  already  mentioned,  chances  are  that  in  a  specific  application  some  of 
the  coefficient  equations  (43)-(48)  are  trivially  satisfied.  In  order  to  obtain  a 
well-posed  problem  it  is  then  necessary  to  eliminate  the  trivial  equations  from 
the  system  and  to  accordingly  reduce  the  matrix  £(w)  as  well  as  the  vectors 
k(w)  and  u(w). 

G.  Formal  solution  and  calculation  of  eigenfrequencies 

The  last  steps  toward  the  solution  of  our  perturbational  problem  are  now 
clearly  prescribed.  Equation  (48)  has  to  be  solved  for  the  unknown  expansion 
coefficients  as  functions  of  u: 

u(«)  =  £-,(«)-Mw)  (50) 

where  j2Tl(w)  is  the  inverse  matrix  to  fi(w).  Application  of  the  inverse  time 
Laplace  transformation  (8b)  to  this  equation  yields  the  time-dependent  coeffi¬ 
cient  vector 

«(*)  =  <r'w‘  =  J  £  e~lut  k(w).  (51) 

i.e.f 

m  -  /  £ >»  (52) 
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- M - 

f  £  e— ‘  £(*,w) 

(53) 

/r(M-j 

(I.-- /fix. -> 

(54) 

7r(M)«  J 

(55) 

The  Anal  solution  functions  are  then  given  by  (52)  and 

£(*.*)-££(*.*)*.(*)  (58) 

* 

(57) 

x 

fUv  <  o,  t)  -  53  t)  Ifi-(v) .  (58) 

a 

As  in  Landau's  infinite-plasma  case,1  some  general  properties  of  these  solu¬ 
tion  functions  can  be  conveniently  retrieved  by  choosing  an  appropriate  contour 
of  integration  in  the  complex  w  plane,  ef.  Eqs.  (52)-  (55)  and  Fig.  2.  Let  the 
“original"  contour  Cx  be  an  infinite  straight  line  with  Imw  »  (Imu)i  **  const., 
where  (Imw)!  must  be  sufficiently  positive  for  C|  to  tie  above  all  singularities 
(poles  and  branch  cuts)  of  the  integrand  vector,  £-1(w)*k(w).  As  is  well 
known  from  function  theory,  the  same  result  u(f)  will  follow  for  any  other 
contour  that  is  obtained  by  deforming  Cx  without  crossing  the  singularities  of 
the  integrand.  Let  us  construct  a  specific  contour  Ct  by  shifting  Cj  down  to 
Imw  -  (Imw)j  <  0  but  excluding  the  singularities  as  shown  in  Fig  2.  The 


solution  (51)  may  then  be  written 


«(*)  —  [-*  Res«^  {ST1  *ke-fe‘)]  +  £  [-»  Res«^/  (l>  1  •  k  e"‘w<)] 

V  V* 


(50) 


where  the  sums  over  v  and  6  respectively  contain  the  contributions  from  the 
poles  w„  and  branch  cuts  Bt  of  £-1(w)  (with  Imw  >  (Imw)j;  RetUr  indicates 
the  residue  at  a  pole  w„),  the  sums  over  t/  and  V  contain  the  contributions 
from  the  poles  w„*  and  branch  cuts  B»»  of  k(w)  (with  Imw  >  (Imw)j),  and 
the  last  integral  represents  the  contributions  from  the  straight  portions  with 
Imw*  (Imw)2. 


As  usual,  we  call  the  system  unstable  if  the  solution  (50)  contains  com¬ 
ponents  growing  in  time  after  the  initial  transient  has  decayed.  Such  con¬ 
tributions  can  arise  from  poles  and  branch-cut  portions  with  Imw  >  0.  An 
nth-order  pole  w*  of  the  r.h.s.  of  Eq.  (50)  gives  rise  to  a  contribution  of  the 
form  P»_i(f)exp(— iuipt),  where  P»_j(t)  is  a  polynomial  of  degree  n  -  1.  In 
particular,  for  the  most  common  case  of  a  first-order  pole  (n  =  1)  we  obtain  a 
purely  exponential  time  behavior. 

The  matrices  £(w)  and  £~l(w)  do  not  depend  on  the  perturbations  but 
are  solely  determined  by  “intrinsic"  equilibrium,  boundary,  and  external-circuit 
properties.  The  poles  w„  of  £-1(w)  are  the  roots  of  the  “characteristic  equa¬ 
tion", 

|£(w„)|«0  (*-...1,2,...)  (60) 

where  (fi(w)(  is  the  determinant  of  Q{ w).  In  accordance  with  common  terminol¬ 
ogy,  the  roots  w„  are  the  eigenfrequencies,  and  the  related  oscillation  patterns 
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'ia  x,  o ,  and  t  are  the  eigenmodes  of  our  system.  Thus,  the  first  term  on  the 
r.h.s.  of  Eq.  (59)  represents  the  eigenmode  contributions  to  the  whole  solution. 

In  a  given  application  it  may  often  be  sufficient  to  study  the  dispersion  and 
stability  behavior,  rather  than  finding  the  detailed  solution  functions  (52)  and 
(56)-(58).  One  then  still  has  to  determine  the  uppermost  poles  and  branch- 
cut  portions  of  both  QT1^)  and  k(w).  However,  in  many  cases  it  will  be 
sufficient  to  consider  the  uppermost  pole,  and  it  may  also  be  possible  to  obtain 
the  determinant  |£(w)|  in  analytical  form,  as  is  the  case  for  the  problem  treated 
in  Part  II.  * 

Due  to  lack  of  sufficient  experience  with  the  present  method,  an  exhaustive 
discussion  of  the  possible  forms  of  behavior  of  these  singularities  cannot  be 
given  as  yet  and  is  therefore  deferred  to  an  appropriate  future  date. 


j 


ra.  DISCUSSION  AND  COMPARISON  WITH  PREVIOUS  LIT¬ 


ERATURE 

A.  Discussion 

The  formalism  developed  in  Sec.  II  is  intended  as  a  first  step  toward  a  com¬ 
prehensive  kinetic  treatment  of  bounded  plasma  systems.  By  ‘‘comprehensive’’ 
we  mean  that  plasma,  boundary,  and  external-circuit  effects  are  to  be  taken 
into  account  simultaneously  and  realistically. 

As  the  first  step  it  represents,  the  present  work  is  still  subject  to  some 
restrictions  which,  however,  may  possibly  be  relaxed  at  a  later  stage,  cf.  be¬ 
low.  The  most  restrictive  assumptions  are  those  concerning  the  linearity  of  the 
problem  (which  means  that  we  restrict  ourselves  to  small  perturbations  about 
a  given  equilibrium  state)  and  the  one-dimensional  system  geometry  involving 
a  uniform  plasma  region  and  thin  sheaths  (whose  relevance  has  been  discussed 
in  Subse.  n.A).  Within  this  framework,  however,  the  admissible  equilibrium, 
initial,  boundary,  and  external-circuit  conditions  are  certainly  suitable  to  model 
a  wide  range  of  realistic  situations.  No  kinetic  treatment  whatsoever  is  known 
to  us  which  includes  these  latter  features  in  comparable  generality,  cf.  Subsec. 

ra.B. 

The  key  result  of  the  present  paper  is  the  system  of  integral  equations 
(37j-(40),  which  is  essentially  equivalent  to  all  the  basic  and  auxiliary  relations 
entering  the  theory  but,  at  the  same  time,  much  more  concise  and  “advanced" 
than  these.  This  is  because  in  deriving  it  we  have  already  performed,  in  general 
terms,  many  of  the  mathematical  manipulations  that  would  be  required  anyway 
when  solving  for  any  non-trivia!  special  case.  We  may  state  that,  within  the 


limits  of  the  model  chosen,  the  treatment  of  the  perturbation  al  problem  has 
been  formalized  to  the  extent  that  seems  possible,  so  that  in  practice  the  “real" 
work  needs  to  set  in  only  at  a  level  where  the  peculiar  details  of  the  special  case 
considered  become  important. 

As  already  mentioned,  we  bear  some  hope  that  it  may  become  possible  at  a 
later  stage  to  relax  the  major  restrictions  still  inherent  to  the  present  work.  This 
hope  is  founded  in  the  fact  that  integral-equation  approaches,  which  tend  to  be 
more  compact  than  the  equivalent  differential-equation  formulations,  seem  to  be 
widely  (if  not  generally)  applicable  to  the  compr  hensive  kinetic  description  of 
bounded  plasma  systems  which  is  attempted  here.  For,  the  sources  of  integral 
terms  spotted  in  the  present  paper,  namely  the  trajectory  integral  (6),  the 
current-conservation  relation  (15)  (with  (IB)),  and  the  boundary  conditions  (28), 
will  also  persist  in  more  general  models.  In  spite  of  the  added  complexity  to  be 
expected,  extensions  of  the  present  work  (e.g,  to  two  dimensions,  non-uniform 
equilibria,  and/or  collisions!  plasma  behavior)  seem  to  be  worthwhile  and  will 
be  attempted.  Specific  examples  involving  the  solution  of  integral  equations  are 
provided  by  Refs.  23-25  (where  the  stability  of  non-uniform  diode  equilibria  was 
studied)  or  Refs.  28  and  27  (where  path-integral  formulations  of  the  linearized 
Boltzmann  equation  were  used  to  solve  collision al  swarm-drift  problems). 

B.  Comparison  with  previous  literature 

In  order  to  locate  the  position  of  the  present  work  within  the  pertinent 
literature,  we  now  present  a  comparison  with  previous  contributions  to  the 
kinetic  theory  of  linear,  longitudinal,  collisionless  plasma  oscillations  in  one¬ 
dimensional  geometry.  An  extremely  large  number  of  relevant  treatments, 
henceforth  referred  to  as  class  A,  has  accumulated  over  the  years.2  Hence, 
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the  list  of  references  to  be  quoted  here  cannot  claim  completeness,  but  we 
nevertheless  believe  it  to  be  representative. 

Let  us  first  note  that  the  overwhelming  majority  of  class-A  treatments  are 
those  concerned  with  Landau’s  infinite-plasma,  initial-value  problem,1'*  cf.  Sec. 

I.  Within  the  framework  of  the  present  method,  this  case  is  the  most  degenerate 
one  and  can  be  readily  retrieved,  e.g.,  from  Eq.  (11)  as  described  there.  What 
then  remains  to  be  considered  is  the  much  smaller  number  of  those  class-A 
treatments  that  actually  involve  boundaries  and  will  henceforth  be  referred  to 
as  class  B. 

Table  II  provides  a  fairly  detailed  comparison  between  17  pertinent  treat¬ 
ments,  namely  Landau's  classical  infinite-plasma  problem  (Ref.  1,  §1),  15  pub¬ 
lished  class-B  treatments,  and  the  preesent  work.  Each  treatment  is  evaluated 
with  respect  to  seven  basic  criteria,  and  the  results  of  the  comparison  may  be 
summarized  as  follows. 

Criterion  1  (model).  Only  nine  references,  including  the  present  work, 
employ  truly  kinetic  theory,  whereas  the  other  eight  are  based  on  the  cold- 
fluid  approach  but  are  still  included  here  because  this  type  of  analysis  is  most 
frequently  encountered  and  easily  reproducible  from  kinetic  theory,  cf.  Part 

II. *  The  fact  that  the  cold-Buid  case  can  be  retrieved  from  the  present  method 
as  a  special  case  is  indicated  by  the  symbol  (+). 

Criterion  2  (number  of  boundaries).  As  the  only  representative  of  the  huge 
number  of  infinite-plasma  treatments*’*  we  have  included  Landau's  classical 
initial-value  analysis  (Ref.l,  §1).  Four  references  are  concerned  with  semi¬ 
infinite  plasmas,  whereas  true  diode  configurations  are  dealt  with  in  twelve 


treatments,  including  the  present  one.  The  latter  coven  the  infinite  and  semi- 
infinite  geometries  as  special  cases. 

Criterion  3  (degree  of  non-uniformity  of  the  d.c.  plasma  equilibrium). 
As  already  stated,  the  present  venion  of  our  method  is  restricted  to  uniform 
plasmas  with  thin  sheaths  and  thus  does  not  include  the  problems  of  Refs.  23, 
24,  25,  31,  and  33  as  special  cases. 

Criterion  4  (description  of  spatial  behavior).  In  five  treatments,  the  spatial 
eigenmode  profiles  are  constructed  by  suitable  superposition  of  infinite-plasma 
modes.  References  13  and  30  start  out  with  general  Fourier-series  expansions, 
whereas  spatial  Laplace  transforms  are  employed  in  Refs,  lg,  19.  20  and  in  the 
present  work.  In  the  remaining  six  references,  the  integral  equations  governing 
the  spatial  eigenmode  structures  ate  directly  derived  from  the  basic  equations. 

Criterion  5  (description  of  time  behavior).  Eleven  of  the  17  references 
considered  only  deal  with  eigenmodes  or  driven  modes  characterized  by  a  time 
dependence  of  the  form  exp(crt).  Reference  33  in  addition  considers  the  effect 
of  an  externally  applied  delta  pulse,  whereas  Laplace  transformations,  covering 
the  full  initial-value  problem,  are  used  in  six  references  including  the  present 
work. 

Criterion  8  (boundary  conditions  on  particles).  None  of  the  previous  treat¬ 
ments  allows  for  boundary  conditions  other  than  constant  emission,  plain  ab¬ 
sorption  (i.e.,  absorption  without  the  release  of  secondary  particles),  and  specu¬ 
lar  reflection.  These  are  all  special  cases  of  the  more  general  relations  (28).  The 
boundary  conditions  assumed  in  Refs.  18  and  20  are  not  immediately  clear. 

Criterion  7  (a.c.  external-circuit  properties).  Ten  previous  treatments  deal 
with  the  external  short  circuit,  and  two  are  concerned  with  the  open  circuit.  No 
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explicit  statement  about  external  circuits  is  made  in  Refs.  1  (§1  and  §2),  13,  18 
19,  20,  and  for  Refs.  13, 18,  and  20  it  would  require  some  detailed  work  to  recon¬ 
struct  the  external-circuit  conditions  that  were  implicitly  assumed.  Interpreting 
Landau’s  infinite-plasma  problem  (Ref.  1,  §1)  in  terms  of  the  diode  picture  as 

~  a 

discussed  after  Eq.(ll),  the  conditions  /  (±oo,  n.f)  =  0  and  V(±oo,t)  =  0  im¬ 
ply  constant  emission  (criterion  8)  and  external  short-circuit  (present  criterion). 
Of  the  treatments  considered,  only  Ref.  31  and  the  present  method  allow  for 
non-trivial  external  circuits  of  some  generality. 

In  summarizing  the  foregoing  comparison  we  may  state  that  the  present 
method  covers  the  problems  treated  in  Refs.  1  (§1  and  §2),  13,  18,  19,  20,  28, 
29,  30,  32,  and  34  as  special  cases.  This  is  not  true  of  Refs.  23,  24,  25,  31, 
and  33  because  these  are  concerned  with  truly  non-uniform  plasma  equilibria. 
Apart  from  this,  all  previous  treatments  are  less  general  than  the  present  one  in 
that  they  start  out  with  more  restricted  specifications  and  thus  cover,  a  priori, 
a  narrower  range  of  the  physical  effects  to  be  expected  in  real  systems.  In  par¬ 
ticular,  it  appears  that  the  effects  of  realistic  boundary  conditions  and  external 
circuits  have  not  been  treated  to  a  significant  extent,  and  that  the  full  pertur- 
bational  problem  for  the  diode  (i.e.,  the  problem  involving  initial  and  external 
perturbations)  has  not  been  theoretically  analyzed  at  all.  We  thus  believe  that 
the  present  method  represents,  in  several  aspects,  a  suitable  starting  point  for 
more  comprehensive,  realistic,  and  systematic  kinetic  descriptions  of  bounded 
plasma  systems  than  have  been  given  so  far. 


IV.  CONCLUSIONS 

We  have  proposed  a  theoretical  method  for  studying  linear  longitudinal 
perturbations  in  one-dimensional  collisionless  plasma  diodes  with  thin  sheaths. 
The  key  result  is  the  set  of  2  4-  2na  integral  equations  (37)-(40)  for  the  2  + 
2n„  time  Laplace  transforms  ;«((>>),  E(x,u),  f*(v  >  0, w),  and  /*(»  <  0, «). 
Upon  expanding  the  functions  involved  in  suitable  basis  sets,  these  equations 
transform  into  the  matrix  equation  (48)  for  the  unknown-  coefficient  vector  u. 
In  principle,  the  full  solution  of  the  perturbational  problem  can  be  found  by 
solving  these  equations,  in  either  representation,  for  the  above  functions  of  u, 
and  by  then  Laplace  inverting  the  latter  into  the  corresponding  time-dependent 
perturbations. 

However,  as  discussed  in  Subsec.  ILC,  the  linear  dynamics  of  the  system 
is  basically  governed  by  the  singularities  of  the  coefficient  determinant  j£(<*>)|, 
so  that  it  may  often  be  sufficient  to  study  these,  rather  than  evaluating  all  the 
details  of  the  full  time-dependent  solution. 

Apart  from  the  geometrical  simplicity  of  the  model,  the  present  method  is 
comprehensive  in  that  it  allows  for  very  general  equilibrium,  initial,  boundary, 
and  external-circuit  conditions  and  thus  combines  a  good  deal  of  the  essential 
elements  present  in  real  systems.  The  literature  survey  of  Subsec.  m.B  suggests 
that  there  exjsts  no  other  kinetic  treatment  which  simultaneously  includes  all 
of  these  aspects  in  comparable  generality. 

Thus,  the  present  method  may  well  represent  a  suitable  starting  point  for 
more  complete,  systematic,  and  detailed  kinetic  studies  of  bounded  plasma 
systems.  A  first  application  is  provided  by  Part  Q  of  the  present  work,*  where 


it  is  demonstrated  that  the  Pierce  instability3*  can  be  drastically  modified  by 
a  non-triviai  external  circuit. 

As  discussed  in  Subsec.  m.A,  the  integral  equations  (37)— ( 40 )  are  equivalent 
to  all  the  basic  and  auxiliary  equations  entering  the  theory,  but  at  the  same  time 
much  more  concise  and  ’’advanced”  than  these.  Due  to  the  relative  geometrical 
simplicity  of  the  one-dimensional  diode  model  considered,  it  has  been  possible 
here  to  perform  a  major  portion  of  the  necessary  mathematical  manipulations 
in  general  terms.  This  means  that  the  whole  approach  has  been  formalized  to 
a  large  extent,  so  that,  in  a  given  application,  the  'real"  work  has  to  set  in 
only  at  a  level  where  the  peculiar  details  of  the  special  case  considered  become 
important. 

According  to  Subsec.  m.A,  it  seems  plausible  that  integral-equation  for¬ 
mulations  of  the  kinetic  perturbational  problem  can  also  be  found  for  bounded 
plasma  systems  that  are  more  general  than  the  one  considered  here.  In  spite  of 
the  added  complexity  to  be  expected,  it  may  well  be  possible  and  worthwhile 
to  extend  the  present  approach  to  some  of  those. 
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APPENDIX:  KNOWN  FUNCTIONS  AND  INTEGRAL  OPERATORS 
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TABLE  I.  Level  scheme  for  establishing  the  coefficient  functions  and  integral 
operators  of  Eqs.  (37H40). 
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TABLE  □.  Comparison  of  some  representative  contributions  to  the  kinetic 
theory  of  linear  longitudinal  plasma  oscillations  in  one-dimensional  geometry. 


Pierce _ 1948 

Rosenbluth  . . .  1963 

present 
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Flgurn  Caption* 

FIG.  1.  Model  geometry,  with  one-minimum  equilibrium  potential  distribu¬ 
tion  as  an  example. 

FIG.  2.  Typical  contours  of  integration  for  the  inverse  time  Laplace  transfor¬ 


mation. 
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